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We consider the Dicke Hamiltonian, a simple quantum-optical model which exhibits a zero- 
temperature quantum phase transition. We present numerical results demonstrating that at this 
transition the system changes from being quasi- integr able to quantum chaotic. By deriving an 
exact solution in the thermodynamic limit we relate this phenomenon to a localisation-delocalisation 
transition in which a macroscopic superposition is generated. We also describe the classical analogues 
of this behaviour. 



At zero temperature, systems of N interacting parti- 
cles can exhibit a quantum phase transition (QPT) as 
a function of a coupling parameter A in the limit that 
N — > oo. How do the precursors of such a transition in- 
fluence quantum chaotic (and non-chaotic) behaviour of 
the same system for finite N? 

One of the most direct indicators of the emergence 
of quantum chaos is the change in energy level spac- 
ing statistics from Poissonian to being described by the 
Gaussian ensembles of Random Matrix Theory. Al- 
though this change-over has been observed in many sys- 
tems 0, 0, 0, , only in a comparatively few, isolated 
cases has the onset of quantum chaos been correlated 
with the presence of a QPT. Important examples in- 
clude spin glass shards, which have recently been used 
in modeling the onset of chaos in quantum computers 
k|, the Lipkin model 6], the interacting boson model 
0, and the three-dimensional Anderson model 0, 
where the change in level statistics occurs at the metal- 
insulator (localisation-delocalisation) transition found in 
disordered electronic systems. 

In this Letter we consider the Dicke Hamiltonian (DH) 
[HI , a quantum-optical model describing the interaction 
of N two- level atoms with a number of bosonic modes. 
We demonstrate that a crossover between Poisson and 
Wigner-Dyson statistics in this model for finite N is in- 
timately connected to a mean-field type superradiance 
QPT. 

The simplicity and generality of the Dicke Hamilto- 
nian have afforded it appeal both for the investigation 
of quantum chaos, and as a model for phase transitions 
at a critical coupling A c induced by the interaction with 
light. The level statistics for finite N have revealed the 
existence of quantum chaos in certain isolated regimes of 
the model |ll|,ll2|. On the other hand, the QPT aspect 
for TV — > oo has been discussed in the context of super- 
radiance 0, Q 5 and recently for exciton condensation 
0. Here, we derive an exact solution for all eigenstates, 
eigenvalues and critical exponents in the thermodynamic 
limit, and show that above the critical point A = A c the 
ground-state wavefunction bifurcates into a macroscopic 
superposition for any N < oo. Our numerical results 
indicate that a localisation-delocalisation transition for 



N — > oo underlies the cross-over between Poissonian and 
Wigner level-spacing distributions for finite N. Further- 
more, we use an exact Holstein-Primakoff transformation 
to derive the classical limit of the model for arbitrary N 
and find a transition at A = A c from regular to chaotic 
trajectories. The latter are delocalised around two fixed 
points in phase space which we conjecture constitute the 
classical analogue of the macroscopic superposition. The 
emergence of this state in the Dicke model is related to 
a conserved parity which becomes spontaneously broken 
in the thermodynamic limit. 

We consider the single-mode Dicke Hamiltonian (DH) 

H = ujqJ z + uoa ] a + — ^= (a) + a) ( J+ + J_) , (1) 

where J z , J± are the usual angular momentum opera- 
tors for a pseudo-spin of length j = N/2, and a, a) are 
the bosonic operators of the field. The atomic level- 
splitting is given by cjo, w is the field frequency, and 
A is the atom-field coupling. Crucially, we have not 
made the rotating-wave approximation (RWA), as this 
would render the model integrable and destroy the phe- 
nomena that we describe here [l2|. There is a con- 
served parity n associated with the DH, which is given 
by n = exp {in [a^a + J z + j] }, such that [H, U] = 0. 
The eigenvalues of n are ±1 and, unless stated, we shall 
work exclusively in the positive parity subspace. 

We begin by discussing the properties of the sys- 
tem in the thermodynamic limit iV, j —> oo. In this 
limit the system becomes integrable for all A, and we 
can derive effective Hamiltonians to describe the sys- 
tem exactly in each of its two phases. We employ a 
procedure similar to Hillery and Mlodinow's analysis of 
the RWA Hamiltonian [l^, and introduce the Holstein- 
Primakoff represent ation of the angul ar mome ntum op- 
erators, J + = tfyj2j - fttfe, J- = y/2j -tfb 6, J, = 
(b^b — j), where b and are bosonic operators 
Making these substitutions allows us to write the DH as a 
two-mode Hamiltonian. Below the phase transition, we 
proceed to the thermodynamic limit by expanding the 
square roots and neglecting terms with powers of j in 
the denominator. This yields the effective Hamiltonian 
= uotfb + uoa^a + A (a f + a) (&t + b) - jco . This 
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FIG. 1: Excitation energies e± of the DH in the thermody- 
namic limit. Inset: scaled ground-state energy, Ec/j, in the 
thermodynamic limit (solid line) and at various finite values 
of j = 1/2, 1, 3/2, 3, 5 (dashed lines). The Hamiltonian is on 
scaled resonance uj = ljq = 1. 
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FIG. 2: Plots of nearest-neighbour distributions P(S) for the 
Dicke Hamiltonian, for different couplings A and pseudo-spin 
j. Also plotted are the universal Poissonian (dots) to Wigner 
(dashes) distributions. 



bi-linear Hamiltonian may then be diagonalised to give 
— £ ^c f _ c _ +£^c+c+ —jouo, where s^± are the exci- 
tation energies of the low-coupling phase, and are given 

by 

2 2 = + u>l ± ^/(^-^ 2 ) 2 + 16A 2 c^ . (2) 

The energy is real only for A < = A c , 

which locates the phase transition. We derive an effec- 
tive Hamiltonian above A c by first displacing each oscil- 
lator mode in the Holstein-Primakoff DH by a quantity 
proportional to ^fj, an d then neglecting terms as above. 
With an appropriate choice of displacements, this pro- 
cess also yields a bi-linear Hamiltonian, which may be 
diagonalised to a form similar to but with different 

vacuum and excitation energies, the latter of which are 
given by 

2A^) 2 = ^ + J\t 

±^( w §A4- w 2A4) 2 +4^ w gAf. (3) 

(2) 

The excitation energy e_ of this second Hamiltonian 
is real only above the phase transition. There are 
two independent choices of displacements, each of which 
yield a different effective Hamiltonian above A c . This is a 
consequence of the fact that at the QPT the n symmetry 
is spontaneously broken. 

In Fig. n we plot as a function of coupling the be- 
haviour of the excitation energies and the ground-state 
energy, which is itself continuous but possesses a discon- 
tinuity in its second derivative at A c . Below the phase 
transition the ground-state is composed of an empty field 
with all the atoms unexcited and hence {J z )g/j = — 1 



and (o)o)q/ j = 0. Above A c , the field is macroscopically 
occupied, (a^a)c/j = 2 (A 4 — A 4 ) / (c^A) 2 , and the atoms 
acquire a macroscopic inversion, (J^g/j = — A 2 /A 2 . At 
the phase transition, the excitation energy S- vanishes as 
|A — \ c \ zv and the characteristic length scale I- = \j -JeZ 
diverges as |A — A c | _zy , with exponents given by v = 1/4, 
z = 2 on resonance. 

To investigate the level statistics of the system, we 
numerically diagonalise the Hamiltonian in the basis 
{|n) \j,m)}, where a^a\n) = n\n), and \j,m) are the 
Dicke states, J z \j,m) = (b^b — j) \j,m) = m\j,m). We 
restrict ourselves to the positive parity subspace by con- 
sidering only states with n + m + j even. We then un- 
fold the resulting energy spectrum to rid it of secular 
variation, form the level spacings S n = E n +i — E n , and 
then construct the nearest-neighbour distribution func- 
tion P (S) . Finally, we normalise the results for com- 
parison with the universal ensembles of Random Matrix 
Theory. In the following, we shall use the term "quasi- 
integrable" to denote systems exhibiting Poissonian level 
statistics, and reserve "integrable" for systems possessing 
exact solutions. 

Figure shows the P (S) distributions obtained for 
the DH at various values of the coupling, and for various 
values of j. At low j the P (S) clearly do not corre- 
spond to any of the universal ensembles. This is most 
keenly observed in the j = 1/2 case (identical to the 
Rabi Hamiltonian 18]), where the spectrum is of "picket- 
fence" character 19], characteristic of one-dimensional 
systems or harmonic oscillators |2£j. For couplings less 
than the critical value, A < A c , we see that as we increase 
j, the P (5) approaches ever closer the Poissonian distri- 
bution, Pp (S) = exp(— S). At and above A c the spec- 
trum is seen to converge onto the Wigner distribution 
Pw (S) = irS/2ex.p (— 7r6' 2 /4), characteristic of quantum 
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FIG. 3: Nearest-neighbour spacing S n = En+i—E n vs. eigen- 
value number n plot for j = 5 with A = 8A C . Horizontal 
crosses: results for the integrable A — > oo Hamiltonian. Inset: 
j = 5 results with A — 4A C and A — 6A C . 
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FIG. 4: The modulus of the ground-state wavefunction 
ip(x,y) of the Dicke Hamiltonian in the abstract x-y repre- 
sentation for finite j = 5, at couplings of A/A c = 0.4, 1.0, 1.2, 
1.4. Black corresponds to Max|-0| and white corresponds to 
zero. 



chaos. This demonstrates that the precursors of the QPT 
in this model lead to a cross-over from quasi-integrable 
to quantum chaotic behaviour at A ~ A c . 

A further transition between integrable and chaotic be- 
haviour is observed in the sequence of level spacings S n . 
The A — > oo limit of the Hamiltonian is integable for 
arbitrary j, having eigen-energies E nm = jn — j^-m 2 , 
where n = 0, 1, 2, . . . and m = — j, . . . , +j. As A is in- 
creased from A c to approach this limit with j fixed, the 
spectrum reverts from Wigner-like to integrable. How- 
ever, it does not follow the usual transition sequence, 
but rather through a sequence illustrated by Fig. [31 The 
spectrum becomes very regular at low energy, where it 
approximates the integrable A — > oo results closely. Out- 
side the regular region the spectrum is well described by 
the Wigner surmise, and the energy-scale over which the 
change between the two regimes occurs is seen to be sur- 
prisingly narrow. As coupling is increased, the size of the 
low-energy integrable window increases, until it eventu- 
ally engulfs the whole spectrum as A — > oo. 

We now proceed to consider the wavefunctions 
of the system by introducing an abstract position- 
momentum representation for each of the boson modes 
via x = -7= (a) + a) ; p x = i^/^(a) — a), and?/ = 

^= (V + 6); p y = iy^f - b). In this representa- 
tion the action of the parity operator n corresponds to 
rotation by tt about the origin. The ground state of the 
system on scaled resonance for j = 5 is plotted for dif- 
ferent couplings in Fig. QJ These wavefunctions were 
obtained by diagonalising the Hamiltonian in the same 
basis as was used to calculate P (5) and representing the 
basis vectors |n) \m,j) by products of harmonic os- 
cillator eigenfunctions. For the non-interacting system 
(A = 0), the wavefunction is a product of two indepen- 



dent Gaussians. As the coupling increases, the two modes 
start mixing, leading to a stretching of the single-peaked 
wavefunction. Around the critical coupling (A w A c ), the 
wavefunction bifurcates to become a double-peaked func- 
tion, and with further increases in coupling the two lobes 
at (±#o, TUo) move away from each other in their respec- 
tive quadrants of the x — y plane. Since xo and yo are of 
the order of y^? f° r large j the ground state evolves into 
a superposition of two macroscopically distinguishable 
parts, which may be considered as a "Schrodinger's cat" . 
The formation of this state constitutes a delocalisation of 
the ground-state wavefunction, which is also observed in 
the excited states. This localisation-delocalisation tran- 
sition is consistent with the transition between Poisson 
and Wigner distributions in the spectrum [2l| . The sup- 
pression of chaos at low j is then seen to be due to the 
fact that for low j only a few excitations are permitted in 
the b- mode. This restricts the extent of the wavefunction 
in the ^-direction, inhibiting delocalisation, and yielding 
the non-generic P (S) seen in Fig. [2j It should be noted 
that an actual spontaneous symmetry-breaking transi- 
tion occurs above A c in the j — > oo limit which removes 
the two lobes so far from one another that they cease to 
overlap and thus can be considered independently. This 
breaks the n symmetry of the model, allowing us to ob- 
tain the earlier exact results by using an effective Hamil- 
tonian for each lobe. These Hamiltonians are identical in 
form, have identical spectra, thus demonstrating that in 
the thermodynamic limit, each energy level in the high- 
coupling phase is doubly degenerate and the macroscopic 
superposition is broken in half. 

Finally, this position representation allows us to study 
the classical analogues of this QPT and the accompa- 
nying onset of chaos in a very natural way. By setting 
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FIG. 5: Typical classical phase space projections (p x = p y — 
0) for the classical Dicke Hamiltonian of Eq. 0) with A/A c = 
1.0, 1.06, 1.2, 1.6, for j = 30.0, uj = ujq = 1. Initial conditions 
were x (0) =2/(0) 1, (0) p y (0) 0. The abrupt 
change to complex motion is observed above A = A c 

the commutators and [y,p y ] to zero, the classical 

Hamiltonian corresponding to Eq. is seen to be 

H c \ = -jVQ + ^ {^W + + V + ^ - U) - (J ) 

+2\y/LJLKxy^l 1 - (c^ 2 + _ . (4) 

Space limitations here only allow us to point out two 
significant features. Firstly, below A c there is only one 
fixed point of the flow, namely x = y = p x = p y = 0. At 
A = A c this situation changes abruptly and two new fixed 
points appear at (x, y) = [±x ch =F2/d) with p x = p y = 0, 
where x c \ and y c \ are approximately equal to the centres 
of the wavefunction lobes £o,2/o f° r l ar g e j- Secondly, 
parametric plots of typical trajectories obtained from Eq. 
(@} (Fig. EJ) demonstrate that the system undergoes a 
rapid change at A = A c from a very simple quasi-periodic 
motion to intricate chaotic behaviour, in agreement with 
the results of the quantum model. Note that the corre- 
spondence between this classical system and the original 
quantum one is significantly greater than previous semi- 
classical treatments and that this holds for any j, 
not just for j — > oo. 

We mention that larger system sizes are required to 
check if there exists a critical level statistics of our model 
at A c 9], and that an examination of the exceptional 
points 6] of this model may yield further insight. Future 
work also includes a study of related models to determine 
the generality of the features described here. 

In summary, we have seen that the P (S) distribution 
of the DH at finite N changes from being Poissonian to 
Wigner at approximately A c , indicating the emergence 
of quantum chaos. The ground- state wavefunction bifur- 
cates at this point, forming a macroscopic superposition. 



The underlying quantum phase transition is reflected in 
the classical model derived here by the appearance of 
two new fixed points at A c , where a transition between 
regular and chaotic trajectories occurs. 
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